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Abstract 

The energy and centrality dependences of charged particle pseudorapidity density in relativistic 
nuclear collisions were studied using a hadron and string cascade model, JPCIAE. Both the rela- 
tivistic p + p experimental data and the PHOBOS and PHENIX Au + Au data at RHIC energy 
could be fairly reproduced within the framework of JPCIAE model and without retuning the model 
parameters. The predictions for Pb + Pb collisions at the LHC energy were also given. We com- 
puted the participant nucleon distributions using different methods. It was found that the number 
of participant nucleons is not a well defined variable both experimentally and theoretically. Thus 
it may be inappropriate to use the charged particle pseudorapidity density per participant pair as 
a function of the number of participant nucleons for distinguishing various theoretical models. A 
discussion for the effect of different definitions in nuclear radius (diffused or sharp) was given. 
PACS numbers: 25.75.Dw, 24.10.Lx, 24.85. +p 
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I. INTRODUCTION 



The main focus of the Relativistic Heavy- Ion Collider (RHIC) at the Brookhaven National 
Laboratory (BNL) is to explore the phase transition related to the quark deconfmement and 
the chiral symmetry restoration. The first available experimental data were the energy de- 
pendence of charged particle pseudorapidity density in central Au+Au collisions at A/i nn =56 
and 130 GeV from the PHOBOS collaboration Q. Soon later, the PHENIX collaboration 
published their data of centrality dependence of the charged particle pseudorapidity density 
in Au + Au collisions at y / i nn =130 GeV 0]. Recently the charged particle pseudorapidity 
density in central Au + Au collisions at y / i nn =200 GeV was also reported by PHOBOS col- 
laboration || . The BRAHMS collaboration published more recently their results of charged 
particle density from Au+Au collisions at i/i nn =130 GeV J|, |J. 




FIG. 1: The energy dependence of the charged particle pseudorapidity density at mid- 
pseudorapidity in relativistic p + p and central A + A collisions. 

It has been predicted that the rare high charged multiplicity in the final state of relativistic 
nucleus-nucleus collisions might indicate the formation of the Quark-Gluon-Plasma (QGP) 
phase in the early stage of collisions ||, |7|, §]. In Ref. ||, the centrality dependence of 
the charged multiplicity has been further proposed to provide information on the relative 
importance of soft versus hard processes in particle production and therefore to provide a 
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means of distinguishing various theoretical models. 

The pQCD calculation with assumption of gluon saturation [JTO , [IT]] (referred to as the 
EKRT model later) was first used to study the centrality dependence of the charged particle 
pseudorapidity density at RHIC. The conventional eikonal approach and the high density 
QCD (referred to as the KN model later) Jl2| were also used to investigate the centrality 
dependence and these two methods surprisingly obtained almost identical centrality depen- 
dence. Recently, authors in |13| reported their study of the same issue from the Dual Parton 
Model. It was found that the experimental observation, the charged particle pseudorapidity 
density per participant pair slightly increasing with the number of participant nucleons, was 
reproduced by 



13], but contradicted the results of [11 
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FIG. 2: The number of participant nucleons < N part > as a function of the percentage of total 
cross section. 



In this paper a hadron and string cascade model, JPCIAE [jlj], was employed to study 
this issue further. Within the framework of this model the experimentally measured energy 
dependence of the charged particle mid-pseudorapidity density per participant pair both 
in relativistic p + p and Au + Au collisions at RHIC energies was reproduced fairly well 
without retuning the model parameters. The predictions for the Pb + Pb collisions at 
LHC energy were also given. Both the PHENIX and the PHOBOS |T5j observations 
that the charged particle mid-pseudorapidity density per participant pair slightly increases 
with the number of participant nucleons could be reproduced fairly well by JPCIAE. In 
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studying centrality dependence the focus was put on the uncertainties in the definitions and 
calculations of the number of participant nucleons. However, this study indicated that it is 
not suitable to use the charged particle mid- pseudorapidity density per participant pair as 
a function of the number of participant nucleons to constrain theoretical models for particle 
production, because the number of participant nucleons is not a well defined physical variable 
both experimentally and theoretically The effect of the different definitions in nuclear 
radius (diffused or sharp) on the number of participant nucleons, the charged particle mid- 
pseudorapidity density, and the charged particle mid-pseudorapidity density per participant 
pair as a function of the number of participant nucleons was discussed. A brief version of 



the part contents of this paper has been published elsewhere [16 
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FIG. 3: The charged particle pseudorapidity density at mid-pseudorapidity in Au + Au collisions 
at i/s nn =l30 GeV as a function of the percentage of total cross section. 



II. MODELS 



The JPCIAE model was developed based on PYTHIA f[? |, which is a well known event 



generator for hadron-hadron collisions. In the JPCIAE model the radial position of a nucleon 
in colliding nucleus A (indicating the atomic number of this nucleus as well) is sampled 
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randomly according to the Woods-Saxon distribution 

Po 



p(r) 



(1) 



1 + exp((r - Ra)/ol) 

where po refers to the normal nuclear density, a ~ 0.54 fm stands for the nuclear diffusion 
edge, and Ra is the nuclear radius of nucleus A. The solid angle of the nucleon is sampled 
uniformly in 4ir. Each nucleon is given a beam momentum in z direction and zero initial 
momenta in x and y directions. The collision time of each colliding pair is calculated under 
the requirement that the least approach distance of the colliding pair along their straight line 
trajectory (mean field potential is not taken into account in JPCIAE) should be smaller than 
\J & tot /ft- Here a to t refers to the total cross section. The nucleon-nucleon collision with the 
least collision time is then selected from the initial collision list to perform the first collision. 
Both the particle list and the collision list are then updated such that the new collision 
list may consist of not only nucleon- nucleon collisions but also collisions between nucleons 
and produced particles and between produced particles themselves. The next collision is 
selected from the new collision list and the processes above are repeated until the collision 
list is empty. 
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FIG. 4: The charged particle pseudorapidity density per participant pair at mid-pseudorapidity 
in Au + Au collisions at y / s nn =130 GeV as a function of the number of participant nucleons, 

< N part >. 

For each collision pair, if its CMS energy is larger than a given cut, we assume that 
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strings are formed after the collision and PYTHIA is used to deal with particle production. 
Otherwise, the collision is treated as a two-body collision |TB|, ITDL EO]. The cut (=4 GeV in the 



program) was chosen by observing that JPCIAE correctly reproduces charged multiplicity 



distributions in AA collisions ||14|| . An important feature of JPCIAE at relativistic energies 
is that QCD parton-parton scatterings are included through PYTHIA |[21|| , which causes 
charged particle yields to increase with collision energy as well as centrality since mini- 
jet production rates increase with energy and number of collisions suffered by participant 
nucleons. It should be noted that the JPCIAE model is not a simple superposition of nucleon- 
nucleon collisions since the rescattering among participant nucleons, spectator nucleons, and 



produced particles is taken into account. We refer to [14] for more about the JPCIAE model. 

Since the number of participant nucleons, N part , plays a crucial role in the presentation 
of PHOBOS or PHENIX data we first make a study on N part . As the direct measurement 
of Np ar t is not available, in the fixed target experiments the number of participant nucleons 
from the projectile nucleus with atomic number A, for instance, is estimated by 

Krt = A*(l-§^), (2) 

beam 

where Ezdc refers to the energy deposited in the Zero Degree Calorimeter (ZDC), dom- 
inated by the energy deposition from the projectile spectator nucleons, and E^ m is the 
kinetic energy of beam . However, in the collider experiments, in order to obtain N part 
one has to relate the measurements to the Monte Carlo simulations. In PHENIX, for in- 
stance, simulations for the response of the Beam-Beam Counter and the ZDC were used to 
calculate N part via the Glauber model 0. In PHOBOS, N part is derived by relating HIJING 



simulations to the signals in the paddle counter |L5|. Therefore, N part here is a model- 
dependent variable. 

Theoretically, the number of participant nucleons in a collision of A + B at impact pa- 
rameter h can be estimated in different ways: 

1. In the geometry method N part {b) reads 

N P art(b) = N p A art (b)+N p B art (b), (3) 

Krt( b ) =PaJ dV9(R A - (x 2 + (6 - yf + z 2 ) l ' 2 )d{R B - (x 2 + y 2 ) 1 ' 2 ), (4) 

KrtiP) =PbJ dV9(R B - {x 2 + y 2 + z 2 ) l ' 2 )d{R A - (x 2 + (b - y) 2 ) 1 ' 2 ), (5) 
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FIG. 5: The number of participant nucleons < N par t > as a function of the percentage of total 
cross section, triangles: diffused nuclear radius, rhombuses: sharp nuclear radius. 

where 8(x)=0 if x <0 and 6(x)=l otherwise, pa (pb) is the nuclear density of the 
projectile (target) nucleus and the nuclear density is normalized to the atomic number. 

2. In the Glauber model, N part (b) is calculated through 

N P art(b) = J d 2 sT A (b- s)[l- exjp(-a in T B (s))) + J d 2 sT B (s)[l - exp(-cr in T A (b - s))) 

(6) 

where <Tj n ~40 mb is the inelastic nucleon-nucleon cross section at RHIC energies and 
Ta (Tb) refers to the nuclear thickness function of nucleus A (B) and is normalized to 

a (b) in. 



3. In the dynamical simulation of A + B collisions, N part (b) can be estimated via counting 
the participant or spectator nucleons and averaging over events simulated at a given 
impact parameter b. However, there are multifarious in simulating models such as: 
FRITIOF @, VENUS [gjj, HIJING [|26|, JPCIAE @, UrQMD j2?|, and AMPT 



23 1, etc. . Not only is the theoretical uncertainty related to the calculation of N part (b) 
in each of the above models large but also the definitions of participant nucleons or 
spectator nucleons are different among each other. We give only a necessary description 
for the following dynamical simulations mentioned in this paper: 
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In FRITIOF [21| and HIJING ||26| , the wounded nucleons, i.e., nucleons which 
suffer at least one inelastic collision, are counted and identified as N part (b). It 
should be pointed out here that in FRITIOF and HIJING the produced particles 
from the string fragmentation do not have rescattering. 

In JPCIAE, we have counted the nucleons involved in at least one inelastic 
nucleon-nucleon collision with string excitation and identified them as N part (b). 



In AMPT [29f|, the spectator nucleons, N spec (b), are counted in the final state 



without rescattering and N part {b) is then calculated through 

N part (b) = (A + B)-N spec (b). (7) 

The spectator nucleons here refer to the nucleons with zero transverse momentum 
and beam energy in the final state of the AMPT simulation. 
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FIG. 6: The charged particle pseudorapidity density at mid-pseudorapidity in Au + Au collisions 
at - v /s nn =130 GeV as a function of the percentage of total cross section, triangles: diffused nuclear 
radius, rhombuses: sharp nuclear radius. 

In the PHENIX or PHOBOS experiment, centrality bin was defined by the cut in the 
particle multiplicity distribution and expressed as the percentage of geometrical (total) cross 
section |2|], g. However, in the theoretical calculation it is more convenient to define the 
centrality by impact parameter b. In order to compare the experimental data of centrality 
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dependence with the theoretical results a relation between g and b is required, which is 
obtained from the definition of the geometrical cross section, that reads 

b = \fg x b max , b max = R A + R B , (8) 

where Ra = 1.12A 1 / 3 + 0.54 fm, for instance, is the radius of nucleus A with diffused 
edge. However, Ra u — 7.5 fm is taken in our calculations since this value (corresponding 
to <j geo =7.2 b) was used in |2| to extract the number of participant nucleons in Au + Au 
collisions. We give in Tab. 1 the mapping between the centrality bin (g bin, selected from 
Tab. 1 in 0) and the b bin according to Eq. (8), and the averaged impact parameter, b, 
over the b bin according to the b 2 law. 



Table 1. The mapping between the g bin 
and the b bin in Au + Au collisions 
(case of diffused nuclear radius) 



centrality bin 


bin of b (fm) 


b (fm) 


below 6* 


below 3.67 


2.45 


10 - 15 


4.74 - 5.81 


5.29 


25 - 30 


7.50 - 8.22 


7.87 


30 - 35 


8.22 - 8.87 


8.55 


35 - 40 


8.87 - 9.49 


9.18 


45 - 50 


10.1 - 10.6 


10.4 



* below 5 for PHENIX. 



As the experimental data were averaged over events in each g bin the theoretical results, 
to be compared with the experimental data, are also averaged over the corresponding b bin. 
We denote the number of participant nucleons after averaging over g or b bin as < N part > 
later. 

III. RESULTS AND DISCUSSIONS 

In Fig. 1(a) the experimental data of charged particle pseudorapidity density per par- 
ticipant pair at mid-pseudorapidity in relativistic p + p (open triangles and rhombuses with 
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error bar) and in central A + A collisions (open circles with error bar for Pb + Pb at SPS 
and full circles with error bar for Au + Au at RHIC) |T], [| were compared with the results 
of the JPCIAE model (full stars for p + p, open squares for Au + Au collisions at RHIC, and 
full squares for Pb + Pb at ^/s nn =17.3 and 5500 GeV). In addition, the results from other 
models were plotted as follows: the dashed and dotted-dash curves are from the HIJING 
model |9j with and without jet quenching, respectively, the solid curve are from the EKRT 



model |I]J . The EKRT results were obtained from || directly, except that the EKRT result 



for Au + Au collisions at y / i nn =5500 GeV was taken from |30|]. Fig. 1 (b) is the same as 
(a) but for A + A collisions only and the vertical coordinate here is the charged particle 
pseudorapidity density itself. One knows from Fig. 1 that both the data of p + p and A + A 
collisions at relativistic energies are also reproduced fairly well by the JPCIAE model. 



A 



3.5 



V 3 

d 



2.5 



Au+Au S =1 30 GeV 




• PHENIX 

▲ JPCIAE(Diffused radius) 
JPCIAE(Sharp radius) 



50 100 150 200 250 300 350 400 



<N port > 



FIG. 7: The charged particle pseudorapidity density per participant pair at mid-pseudorapidity in 
Au + Au collisions at i/s nn =130 GeV as a function of < N par t >■ The circles with shaded area 
of systematic errors: PHENIX data, triangles: JPCIAE (diffused radius), rhombuses: JPCIAE 
(sharp radius). 



In Fig. 2 the < N part > extracted by PHENIX and PHOBOS |5f from Au + Au 
collisions at y / i nn =130 GeV were compared with the model computations. The horizontal 
axis in Fig. 2 is the percentage of geometrical cross section, g, and each g bin is represented 
by its middle point for convenience in plotting (the same for Fig. 3, 5, and 6). The solid 
and open circles with error bar in Fig. 2 were PHENIX and PHOBOS results, respectively. 
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< Np ar t > from the geometry method, open stars, were from the geometric N part (b) (Eq. 3 
- 5) after averaging over b sampled randomly in each b bin due to the b 2 law. Similarly, the 
FRITIOF, JPCIAE, and AMPT [g^] results, full stars, full triangle-ups, and full triangle- 
downs in Fig. 2, were obtained averaging over events simulated for b sampled randomly in 
each b bin due to the b 2 law, respectively. In the EKRT model [11|, the curves of N part (b) 
vs. b/R,A were given and N part is calculated by the Glauber method with Ra=1A2A 1 ^ 3 - 
0.86A- 1 / 3 . The open triangle-ups in Fig. 2 are the EKRT model results taken from the 
\/s nn =i30 GeV one of those curves according to the b in Tab. 1. Originally, HIJING N part 



| 3Tfl were calculated for individual 6, the HIJING points, open squares, in Fig. 2 were plotted 
after relating b to g according to Eq. (8). The full squares in Fig. 2 were the results of the 
KN model taken from Tab. 1 in |12[ under the centrality bins of - 6, 10 - 20, 20 - 30, 30 
- 40, and 40 - 50 %, respectively. In the KN model, < N part > was the result of N part (b) 
from the Glauber model after averaging over centrality bin. In this approach, the particle 
multiplicity and impact parameter were related by a Gaussian distribution with parameters 
fixed via fitting the PHOBOS charged multiplicity distribution. As proved in |3"2" |, such kind 
of average is approximately equivalent to the average method based on Eq.(8). From Fig. 
2 one knows that the discrepancies among the PHENIX or PHOBOS and the model results 
are visible. 

The charged particle pseudorapidity density at mid-pseudorapidity in Au + Au collisions 
at -^^=130 GeV as a function of the percentage of geometrical cross section was given 
in Fig. 3. In this figure, the full and open circles with error bar are the PHENIX and 



PHOBOS data, respectively. The full and open triangles, respectively, are the JPCIAE 
results with and without rescattering. One knows from Fig. 3 that the rescattering only leads 
to a few percent increase in the charged multiplicity although rescattering might enhance 
yields of strangeness, S~ + S~ for instance, by a couple of times. 

In Fig. 4 (a) we compared the PHENIX data of charged particle mid-pseudorapidity 
density per participant pair (full circles with shaded area of systematic errors) with 
the results of the JPCIAE model (full triangles) and the results of other models (obtained 
from H directly): HIJING (the dotted curve), the KN model (the solid curve), and EKRT 
(the dashed curve). One sees that except EKRT, three other models predict an increase of 
(dN c h/drj | r,=o)/(0. 5 < Np ar t >) as a function of < N part > though the theoretical results seem 
to underestimate the PHENIX data. Such an increase can be understood in JPCIAE as a 
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result of increasing hard parton scatterings per participant nucleon. Fig. 4 (b) compared the 
PHENIX data to the results of single dN ch /dr}\ v=0 from JPCIAE normalized by the < N part > 
from different models (taken from the corresponding curve in Fig. 2 at the middle point of 
g bins for KN, HIJING, and EKRT models): full squares, < N part > from the KN model, 
open squares from HIJING, open triangles from EKRT, full stars from FRITIOF, and full 
triangles from JPCIAE. One sees from Fig. 4 (b) that starting from the charged particle 
mid-pseudorapidity density obtained in JPCIAE, but using < N part > from different models, 
the corresponding results of (dN ch /drj\ v=0 )/ (0.5 < N part >) are different visibly among 
each other, in peripheral collisions especially. Therefore it might be inappropriate using 
(dN ch / dr)\ v=0 ) / (0.5 < N part >) as a function of < N part > to distinguish various theoretical 
models for particle production since < N part > is not a well defined physical variable. 

If one uses the sharp nuclear radius, Ra = 1.12A 1 / 3 fm (Ra u =Q-50 fm), instead of diffused 
one used above, the corresponding mapping between the centrality bin (g bin) and the b bin 
from Eq. (8) is given in Tab. 2 . Fig. 5, 6, and 7 give, respectively, the comparisons between 
the JPCIAE results of diffused nuclear radius and the JPCIAE results of sharp nuclear radius 
in the < N part > as a function of centrality, the dN c h/ dri\, n= Q as a function of centrality, and 
the (dNch/dril^o) / (0.5 < N part >) as a function of < N part >. One knows from Fig. 5 
and 6 that there are observable discrepancies between the results of diffused nuclear radius 
and the results of sharp nuclear radius , in peripheral collisions especially. However, in the 



Table 2. The mapping between the g bin 
and the b bin in Au + Au collisions 
(case of sharp nuclear radius) 



centrality bin 


bin of b (fm) 


b (fm) 


below 6* 


below 3.20 


2.13 


10 - 15 


4.11 - 5.03 


4.59 


25 - 30 


6.50 - 7.12 


6.82 


30 - 35 


7.12 - 7.69 


7.41 


35 - 40 


7.69 - 8.22 


7.96 


45 - 50 


8.72 - 9.19 


8.86 



* below 5 for PHENIX. 



12 



(dN c h/drj\ v= o)/ '(0.5 < Np art >) vs. < Np ar t > plot both results are close to each other as 
shown in Fig. 7 . That indicates again that it is hard using (dN c h/ dr]\ v =o) I '(0-5 < N part >) 
as a function of < N part > to distinguish various theoretical models for particle production 
since N part is not a well defined variable. As one sees from Fig. 5 and 6 that for most 
central collision both diffused and sharp nuclear radii work well in mapping centrality bin 
and b bin therefore we do not redraw a figure as Fig. 1 for sharp nuclear radius. However, it 
should be pointed out that using sharp nuclear radius in mapping centrality bin and b bin is 
inconsistent with that the nuclear diffusion edge is introduced in Woods-Saxon distribution 
initiating the nucleons in nucleus. 

IV. CONCLUSIONS 

In summary, we used the hadron and string cascade model, JPCIAE, to investigate 
the energy and centrality dependences of charged particle pseudorapidity density at mid- 
pseudorapidity in relativistic p + p and A + A collisions. Both the relativistic p + p experi- 
mental data and the PHOBOS and PHENIX data of Au + Au collisions at RHIC energies 
could be reproduced fairly well within the framework of the JPCIAE model without retuning 
any parameters. The JPCIAE model predictions for Pb + Pb collisions at the LHC energy 
were also given. Both the PHENIX |2| and the PHOBOS Jl5| observations that the charged 
particle mid-pseudorapidity density per participant pair slightly increases with the number 
of participant nucleons could be reproduced fairly well by JPCIAE. This study shows that 
since < N part > is not a well defined physical variable both experimentally and theoreti- 
cally it may be hard to use charged particle pseudorapidity density per participant pair at 
mid-pseudorapidity as a function of < N part > to distinguish various theoretical models for 
particle production. A discussion for the effects of the different definitions in nuclear radius 
(diffused or sharp) is given and it is indicatd that using diffused nuclear radius might be 
better. 
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